Correlation length of the two-dimensional Ising spin glass with Gaussian interactions 
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We study the correlation length of the two-dimensional Ising spin glass with a Gaussian distri- 
bution of interactions, using an efficient Monte Carlo algorithm proposed by Houdayer, that allows 
larger sizes and lower temperatures to be studied than was possible before. We find that the "effec- 
tive" value of the bulk correlation length exponent v increases as the temperature is lowered, and, at 
low temperatures, apparently approaches —1/0, where ~ —0.29 is the stiffness exponent obtained 
at zero temperature. This means scaling is satisfied and earlier results at higher temperatures that 
find a smaller value for v are affected by corrections to scaling. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 
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I. INTRODUCTION 

There are two main theories to describe the spin-glass 
state: the droplet theorji2i2ii and the replica symmetry 
breaking theorj&SiSi (RSB) of Parisi. According to the 
droplet picture the lowest energy excitation, or "droplet," 
of linear size I containing a given site has a characteristic 
energy of order l e where 9 > is a "stiffness" exponent. 
Droplets are expected to be compact but with a surface 
that has a nontrivial fractal dimension d s , less than the 
space dimension d. It is further assumed that the same 
exponent 9 describes both droplet and "domain-wall" ex- 
citations. In the alternative RSB scenario, the energy of 
droplets containing a finite fraction of the system does 
not increase with increasing system size. Furthermore, 
in RSB the surface of the large-scale, low-energy excita- 
tions are expected to fill space and so have the fractal 
dimension d s = d. 

There have been many numerical studies in 
£k reQ 9 4 io i ii 4 i2 i i3 anc j fouri^i^ dimensions that at- 
tempt to determine which of these scenarios, or possibly 
something els o 16 ' 17 , is correct. These calculations 
are quite limited in the range of sizes that can be 
studied, although recentlyi&i 9 " a larger range has been 
studied for a one-dimensional model with power law 
interactions. Another case where a large range of sizes 
can be studied is the two-dimensional spin glass with 
short-range interactions, for which there is no spin-glass 
order at finite temperature, corresponding to 9 < 0. 
It is desirable to understand fully the two-dimensional 
spin glass, including the nature of corrections to scaling, 
since this may help in the interpretation of numerical 
data in higher dimensions. 

However, even in d = 2, the situation is not completely 
clearcut. Zero-temperature calculations of the energy of 
a domain wall 9 i 20 i 21 i 22 i 23 i 24 consistently give 9 ~ —0.29, 
with Ref. |M for example, quoting 9 = -0.287 ± 0.004. 
However, some calculations of droplet energies find 2 ^ 2 ^ 
9 ~ —0.47, while others^ 7 , find results consistent with the 
domain- wall value. These discrepancies presumably arise 
because some of the results are affected by corrections to 



scaling, and, as discussed by Hartmann and Moore^ the 
domain-wall value, 9 ~ —0.29, seems to be the correct 
asymptotic result. 

Although the discrepancy between the estimates for 
9 from the zero-temperature calculations seems now to 
be resolved^ there is still an apparent conflict between 
9 and the value of the correlation length exponent v = 
2.0 ± 0.2 obtained from finite temperature Monte Carlo 
simulations * 29 Since the spin-glass transition occurs at 
T = in two dimensions, the correlation length £ diverges 

a g J 1 > aS £ ' T ~ V livmrJmn 9.10 

related to 9 by 



T v . According to scaling, 9, v is 



(1) 



which gives v ~ 3.5, significantly different from the result 
v = 2.0 ± 0.2 reported in Ref. We investigate this 
discrepancy here by performing Monte Carlo simulations 
on the Ising spin glass with Gaussian interactions in two 
dimensions at larger sizes and lower temperatures than 
in Ref. |H 



II. MODEL AND OBSERVABLES 



The Hamiltonian is 



H 



ijSiSj, 



(2) 



where the sum is over nearest neighbor pairs of sites on 
a square lattice, the Si are Ising spins taking values ±1, 
and the Jij are Gaussian variables with zero mean and 
standard deviation unity. The square lattice contains 
N = L x L sites with periodic boundary conditions. We 
use the Monte Carlo algorithm of Houdayer— which com- 
bines single spin flip dynamics, parallel tempering^ and 
a type of cluster move, which is significantly more effi- 
cient than parallel tempering in two dimensions for large 
system sizes (L > 24). Tests for equilibration are done 
as in Ref. l32l the parameters used in the simulations are 
shown in Table [I] 



TABLE I: Parameters of the simulations. N samp is the num- 
ber of samples, N BVCCp is the total number of Monte Carlo 
sweeps for each of the 2Nt replicas for a single sample, T m i n 
is the lowest temperature simulated, and Nt is the number 
of temperatures used in the parallel tempering method. Note 
that for L < 16 standard parallel tempering Monte Carlo was 
used, whereas for L > 32 the cluster method by Houdayer 



was 


applied. 
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The main focus of our study is the correlation length 
of the finite system 33 ^ 35 ^ 36 £l , denned by 



CI 



1 



2sin(|k min |/2) 



Xsg(O) 

XSG(kmin) 



"I 1/2 



(3) 



where k m ; n = (2ir/L,0, 0) is the smallest nonzero wave 
vector, and 



i,3 



(4) 



is the wave vector-dependent spin-glass susceptibility. In 
Eq. Q the angular brackets (• ■ •) denote a thermal aver- 
age while the rectangular brackets [• • -] av denote an av- 
erage over the disorder. 

Since the ratio £,l/L is dimensionless, it satisfies the 
finite size scaling f rm22i2&i2& 



(5) 



assuming a zero-temperature transition, where X is a 
scaling function and X/v = —8. For L X I V T 1, £l is 
equal to the bulk (i.e., infinite system size) correlation 
length £00, and so 



& = 6» ~ T~ v 

implying that X(x) ~ 
limit, x <C 1, we expect X(x) 
mate A below. 



(L > Coo) , 



(6) 



for x ^> 1. In the opposite 
where we will esti- 



III. RESULTS 

Figure ^ shows data for We see that the data are 
independent of system size at high T, showing that the 
bulk behavior has been obtained, and the data "peels 
off" from this "bulk curve" at temperatures that become 
successively lower for larger sizes. The bulk results are 
curved on this log-log plot showing that the asymptotic 
power law behavior in Eq. © has not yet been reached. 
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FIG. 1: A log-log plot of the finite size correlation length £l 
for different system sizes and temperatures. The data labeled 
64 — > 00 and 128 —* 00 come from an extrapolation to the 
thermodynamic limit. The slope of the extrapolated data 
gives —v — —3.45, which is consistent with 9 = —1/v ~ 
—0.287 found in zero-temperature domain-wall calculations 
(Ref.|H. 



Rather, the slope of the curve is an "effective exponent," 
i/ e ff, which varies with T. 

In order to determine the asymptotic value of v we 
obtain values for the bulk correlation length at lower T 
following the method used by Kim^ The finite-size scal- 
ing expression, Eq. (J5J, can be written as 



_ , ( Coo 

L ~ \L 
which can be inverted to give 



5 ~L 



(7) 



(8) 



where g(x) = f^ 1 {x). Clearly g(x) — x for x — > 0. We 
determine g(x) by fitting to data in the range 0.45 < 
T < 1.05 where we have data for the correlation length 
in both the bulk and finite-size regimes. We consider 
sizes 16 < L < 128 for this determination, from which 
we obtain g(x) in the range < x < 0.45. Using g{x) 
in this range we then determine foo from Eq. (jSJ) using 
data for L = 64 in the range 0.285 <T< 0.482 and for 
L = 128 in the range 0.24 <T < 0.38. Note, that we do 
not perform any direct extrapolation in this analysis; the 
function g(x) is determined by fitting and is then used 
to get Coo at somewhat lower temperature using only the 
range of x where it was fitted. 
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FIG. 2: Scaling plot of the data for the correlation length 
according to Eq. with 6* = — 1/v — —0.29. The dashed 
line, which fits the data at low T, has slope —1/2. 



FIG. 3: Scaling plot of the data for the correlation length £l 
according to Eq. Jjj), with 6 — —l/v = —0.45. 



The resulting values of ^ from the L — 64 and 128 
data are consistent with each other where they overlap, 
and are shown in Fig. The extrapolated data fit a 
slope of — v = —3.45 which corresponds to 9 = — 1/v = 
—0.29 in good agreement with domain-wall results^ At 
higher temperature, the slope of the bulk data in Fig.QJis 
smaller in magnitude, so we believe that Liang's result, 
v = 2.0 ± 0.2, obtained in the region around T = 1.0, 
is only an effective exponent. Neither our results nor 
the domain-wall results appear to be consistent with an 
exponential divergence of the correlation length as T — > 
0. 

Assuming that the asymptotic value of v is indeed ss 
3.45 we estimate from Fig. ^that one needs to be at or 
below a temperature of around 0.35, where £oo ~ 40, to 
see the bulk asymptotic behavior. 

Further insight is obtained by scaling the data ac- 
cording to Eq. JSJ. Figure El shows a scaling plot with 
6 = —0.29, the value expected from zero-temperature 
domain-wall calculations and our extrapolated data in 
Fig. ^ We see the data in Fig. [3 scale well at low T 
but not well at high T. However, the latter point is not 
surprising in view of Fig. ^ where we see that asymp- 
totic bulk power-law behavior has not yet been reached. 
Note, though, that the data for the two largest sizes in 
Fig. |2 L = 64 and 128, do almost collapse, suggesting 
that 9 = —1/v = —0.29 will work in the bulk region, 
L ^> £, for large enough sizes and low enough tempera- 
tures, as we also inferred from Fig. ^ 

The dashed line in Fig. [21 has slope —1/2, implying 
from Eq. © that X(x) ~ x~ x for x < 1 with A ~ 1/2. 



Hence, we have 

& ~ T -i/2 L i-i/(2u) (9) 

The T -1 / 2 dependence can be understood from Eq. © 
since Xsg(0) = L 2 at T = (because the ground state 
is unique) , while the fluctuations at nonzero k are frozen 
out at T = 0. It is plausible XSG(k m in) k T at low T 
from equipartition, and this leads to a T~ x l 2 dependence 
for £ L . 

Figure[3]shows a scaling plot for 9 — —0.45, which gives 
the best data collapse in the high-T region. This value is 
compatible with v — 2.0±0.2 found in Ref.|23. Note, that 
the data do not collapse at all in the low-T region and the 
collapse becomes worse for larger sizes. Figure shows, 
again, that an effective value of v w 2 will fit the data 
over a range of intermediate temperatures, while in the 
low-T asymptotic region one has 9 = — 1/v ~ 0.29. We 
have found that the spin glass susceptibility shows similar 
behavior. Independent recent results for the spin glass 
susceptibility and other quantities^ also find evidence 
that 9 ~ —0.29 for large sizes. 



IV. CONCLUSIONS 

To conclude, we have shown that the result v = 
2.0 ±0.2 in Ref. |2^is only an effective exponent, and the 
true value for v is larger. The data are consistent with 
scaling holding asymptotically for v = — 1/6 ~ 3.45. We 
have strengthened our argument for this conclusion by 
the extrapolation to L = oo shown in Fig. ^ Of course it 



4 



would also be desirable to extend the data to larger sizes, 
which may be possible in the near future by fine tuning 
the algorithm. 
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